Self-optimization of optical confinement in ultraviolet photonic crystal slab laser 
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We studied numerically and experimentally the effects of structural disorder on the performance 
of ultraviolet photonic crystal slab lasers. Optical gain selectively amplifies the high-quality modes 
of the passive system. For these modes, the in-plane and out-of-plane leakage rates may be au- 
tomatically balanced in the presence of disorder. The spontaneous optimization of in-plane and 
out-of-plane confinement of light in a photonic crystal slab may lead to a reduction of the lasing 
threshold. 



Over the past few years, tremendous progress has been 
made in design and fabrication of photonic crystal slab 
(PhCS) lasers that operate in infra-red (IR) spectrum 
range Q, S H S H S 0] ■ To realizc an ultraviolet (UV) 
PhCS laser, the feature size has to be reduced roughly by 
a factor of four 0,0 compared to IR PhCS. Fabrication 
of such fine structures inevitably generates random devi- 
ations from the perfectly ordered structures. In passive 
PhC devices such uncontrollable disorder has detrimental 
effect as it contributes to optical losses and limits light 
propagation length. However, it is not clear how the dis- 
order would affect the performance of an active device, 
e.g., a PhCS laser. In this letter we illustrate, through 
numerical simulation and experiment, that structural dis- 
order in a PhC laser may not be as detrimental as in a 
PhC waveguide. Under some circumstances structural 
disorder enables self-optimization of optical confinement 
in a PhCS, leading to a reduction of lasing threshold. 

A PhCS utilizes index guiding to confine light to the 
plane of the slab ^3 • In-plane confinement is realized ei- 
ther via a defect state located inside a photonic bandgap 
(PBG) 0, 0, SJIIlOT a bandedge state with vanishing 
group velocity [& Q> EH 113 Light may escape from 
the PhCS vertically through the top/bottom interfaces 
into air/substrate or laterally via the edge of the peri- 
odic pattern into air or unpatterned part of the slab. 
The vertical leakage rate is characterized by the out-of- 
plane energy loss per optical cycle QJ 1 , and the lat- 
eral by Q7\ ■ A defect state spatially localized in the 
vicinity of an intentionally-introduced structural defect 
typically has large leakage in the vertical direction, i.e., 
Q± ^ QV\ ■ For a bandedge state, the lateral leakage 
usually dominates over the vertical one, 



The total loss is described by Q t 



Low 



Hot — Q\_ + Q\\ ■ 
lasing threshold demands maximization of Qtot, which is 
hindered by Q± for a defect state and Q\\ for a band- 
edge state. Several designs aim at optimization of PhCS 
lasers by balancing Q \ and Q\\ via "gentle localization'' 
0, e.g., phase-slipjj 0| , double-heterostructure 0. 

Recently we realized the first UV PhCS laser 8]. ZnO 
films were grown on sapphire substrates by plasma en- 



hanced MOCVD |l4j. Hexagonal arrays of cylindrical 
air voids were patterned in the ZnO films by focused 
ion beam (FIB) etching technique. The lattice constant 
a ~ 120 nm, the radius of air cylinders R ~ 30 nm. 
Post thermal annealing was employed to remove the FIB 
damage. Single-mode lasing at room temperature was 
realized with optical pumping. The scanning electron 
micrograph (SEM) of a ZnO PhCS is shown in Fig. 1(a). 
Despite the long-range periodicity exhibited in the in- 
set of Fig. 1(a), Fig. 1(b) reveals the deviation of the 
fabricated pattern from the ideal honeycomb structure. 
Such "crescent" deviation 0] caused optical scattering 
on the length scale of a few lattice constants. It was ex- 
pected to enhance radiative leakage of a PhCS laser based 
on either defect state or bandedge mode. Moreover, the 
propagation loss in a passive PhCS caused by random [lij 
scattering was predicted to increase dramatically near a 
photonic bandedge [l^ . where the bandedge- type PhCS 
laser operates. Despite of these pessimistic expectations 
based on passive systems, we show that the performance 
of a PhCS laser may be less susceptible to the detrimen- 
tal effects of structural disorder. This is because optical 
gain predominantly amplifies the mode with the high- 
est quality factor Qtot- For the highest-Qt ot mode, the 
vertical and lateral leakage rates may be automatically 
balanced in the presence of disorder. This implies that 
an appropriate amount of structural disorder could lead 
to spontaneous optimization of in-plane and out-of-plane 
confinement of light in a PhCS. 

To investigate how the disorder affects the rates of ver- 
tical and lateral leakages of light from a PhCS, we con- 
sider a system schematically depicted in the inset of Fig. 
2(a). A dielectric slab of thickness 180 nm and refractive 
index n = 2.35 is sandwiched between air and substrate 
(n su b — 1-78). Within the slab, N infinitely- long grooves 
run parallel to y-axis. The width of a groove is 22 nm, 
the lattice constant of the disorderless structure is 100 
nm. We consider light propagating in x-z plane, with 
the electric field along y-axis. Such system is 2D, that 
allows numerical simulation of large statistical ensembles 
of random systems. Despite the simplification, the sys- 
tem in Fig. 2(a) retains the property essential for our 
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FIG. 1: (color online) (a) Top-view SEM of a ZnO PhCS. 
The inset shows structural Fourier transform from the digi- 
tized SEM. Long-range periodicity is reflected in the six max- 
ima at the positions corresponding to the perfect lattice, (b) 
Difference between the digitized SEM of real sample and the 
perfect honeycomb lattice reveals the structural disorder. 



study of PhCS laser: the possibility of vertical (along 
z-axis) and lateral (along x-axis) radiative escape. Us- 
ing the finite-difference time-domain (FDTD) method, 
we find the mode of the passive system that has the high- 
est Qtot A Gaussian pulse was launched at all spa- 
tial points in the slab and the energy is allowed to leak 
out radiatively. Simulation area is terminated by uni- 
axially perfectly matched absorbing layer that absorbs 
all outgoing waves. The pulse excites all modes within 
30 nm wavelength range around 400 nm. After the ini- 
tial multi-mode decay the field distribution is stabilized 
and the longest-lived mode can be seen. This is further 
confirmed by observing a mono-exponential decay of the 
total energy 0, 0] stored in the system that allows de- 
termination of Qtot- By integrating Poynting vector over 
the corresponding interfaces [lj, we obtained the outgoing 
flux in the vertical and horizontal directions, and Q± and 

'I 1 + Q^ 1 relation was 
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Q||. In our simulation, Q~ t 
satisfied numerically to within 2%. 

Fourier transform of the spatial profile of electric field 
at the interface between the slab and substrate gives 



the mode's distribution in fc|| (in-plane component of the 
wavevector) space. In a perfectly periodic structure, the 
bandedge mode has the highest-<5 tot . It is spatially ex- 
tended in x, thus has a narrow distribution in fc|| [thick 
dashed curve in Fig. 2(a)]. Next we intentionally create 
a defect by increasing the spacing between two neighbor- 
ing grooves at the center of the pattern to 150 nm. The 
highest-Q tot mode is localized around this artificial defect 
with a localization length of 140 nm. Strong localization 
in x results in a broad distribution in fc|| [thin dashed 
curve in Fig. 2(a)], with the maximum lying closer to 
the edge of substrate light-cone [dash-dotted vertical line 
in Fig. 2(a)]. Its Qtot is limited by Q±, which is about 
three times smaller than the corresponding Qi in a sys- 
tem of N = 24. In contrast, the bandedge mode is con- 
centrated well beyond the light-cone in fc||-space, thus its 
Q± is much higher. However, its spatial extension makes 
the lateral leakage larger, hence its Q to t is limited by Q\\. 

To simulate the position disorder of air cylinders in 
real structure [Fig. 1(b)], random variation of groove po- 
sition x n is introduced. We choose Aa: n randomly from a 
uniform distribution with the standard deviation Sx =5, 
10, 15 nm. 5x characterizes the "strength" of disorder. 
As the disorder is introduced, the highest-Q tot state dif- 
fers from realization to realization, and the correspon- 
dent <3||, Q± as well as the frequency vary. We study 
statistical distributions of these parameters and their de- 
pendences on disorder strength Sx and system size N. 

In small systems (N = 12 and 24) with an artificial 
defect and weak disorder (Sx = 5 nm), the highest-Q tot 
modes always concentrate around the defect at the cen- 
ter of the pattern. These modes become more spatially 
extended than those without disorder. Therefore, their 
k\ | -distribution is narrowed and hi component within 
the light-cone is significantly reduced [Fig. 2(a)]. This 
reduction leads to a decrease in the vertical leakage, 
thus, an increase in Q±. Meanwhile, Q|| starts increas- 
ing as the mode gets less localized in real space. The 
ensemble-averaged (Q\\/Qj_), shown in Fig. 2(b), de- 
creases monotonously to unity with increase of disorder 
strength. Therefore, disorder removes the imbalance be- 
tween vertical and lateral leakages of a single defect state, 
making ~ (Qj_)- As a result, the ensemble-averaged 
quality factor (Qtot) is slightly higher than that without 
disorder. In a larger system or with stronger disorder, 
the highest-Qtot mode is no longer pinned at the artifi- 
cial defect. Instead, it can explore the entire pattern to 
find the optimum configuration for the best vertical and 
lateral confinement. This leads to a further increase of 

(Qtot)- 

With the introduction of disorder, bandedge mode be- 
comes less extended. As its "tail" moves away from the 
boundaries of the pattern, the lateral leakage decreases, 
thus Q\\ increases. Meanwhile, the distribution in k\\- 
space is broadened and shifted closer to the light-cone 
edge [Fig. 2(a)]. The increase in vertical leakage results 
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FIG. 2: (color online) (a) fcjj-distributions of the highest- 
Qtot modes at one pixel beneath the slab/substrate interface. 
Thin/thick dashed curve represents the mode found in the dis- 
orderless system (N = 24) with/without an artificial defect. 
The corresponding solid curves are representative examples 
of the highest-Qtot modes in these systems with position dis- 
order (Sx = 10 nm). The vertical line marks the substrate 
light-cone boundary. The inset is a schematic sketch of the 
simulated structure. Arrows show the directions of in-plane 
and out-of-plane radiative leakages, (b) Squares/circles rep- 
resent (Q\\/Q±), averaged over 300 random realizations of 
TV = 24 system with/without the artificial defect, versus dis- 
order strength Sx. 



in a decrease of Q±. The ensemble-averaged (Q\\/Q±), 
shown in Fig. 2(b), rises continuously to unity with in- 
creasing disorder strength. Again, disorder balances the 
vertical and lateral leakages of the bandedge mode, as 
it does to the defect state. However, for a bandedge 
mode the increase in (Qm) is not as large as the de- 
crease in (Q±), thus (Qtot) is slighter lower than that 
without disorder. Nevertheless, as the pattern size N in- 
creases, the total leakage rate decreases monotonically: 
(Qtot) a N~ a [Fig. 3(a)]. The exponent a decreases 
from 2.3 at 5x = 5 nm to 1.9 at dx — 15 nm. Even with 
the largest disorder we simulated (Sx = 15 nm), no sat- 
uration of (QYot) with N is observed up to N = 48. This 
behavior differs fundamentally from that of a photonic 
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FIG. 3: (color online) (a) (1/Qtot) for the highest-Q tot modes 
found numerically in PhCS depicted in the inset of Fig. 2(a). 
The average is taken over an ensemble of 300 random real- 
izations. The squares and triangles represent the results for 
Sx — 5 nm and 15 nm, respectively. Circles correspond to the 
disorderless system. Dotted, solid and dashed lines are N~ a 
fits with q = 2.7, 2.3 and 1.9 respectively, (b) Ensemble- 
averaged incident pump pulse energy at the lasing thresh- 
old, measured in the samples of a — 115 nm (squares) and 
a — 118 nm (triangles), versus the pattern size. The data are 
fitted with the power law dependence as in (a) with a — 2.5 
(solid line) and 1.7 (dashed line), (c) Intensity profile of the 
highest-quality mode found numerically in a 2D disordered 
hexagonal array of air holes in the dielectric with refractive 
index n = 1.55. The pattern size is 5 x 5 /xm, lattice constant 
a — 186 nm, and the hole radius R = 65 nm. The wavelength 
of the mode is 393 nm. It is located at the dielectric bandedge, 
and the electric field is concentrated inside the dielectric, (d) 
Measured intensity distribution of the lasing mode in a ZnO 
PhCS with a — 115 nm and R — 0.25a. Dashed box outlines 
the boundary of the 8x8 /im pattern. 



crystal waveguide, where optical loss increases exponen- 
tially with its length. In contrast, a disordered PhCS 
laser benefits from an increase of the pattern size, sim- 
ply because a larger system provides a bigger pool of 
modes from which the highest-Qtot mode can be selected. 
This effect should be more pronounced in PhCS micro- 
lasers with 2D periodicity [Fig. 1(a)], due to the larger 
phase-space compared to the numerically simulated sys- 
tems with ID periodicity. 

Experimentally, we fabricated ZnO PhCS of dimen- 
sions 4x4 jttrn, 6x6 /urn, and 8x8 ^m [Fig. 1(a)]. Since 
the complete photonic bandgap in ZnO PhCS without 
"undercut" was quite narrow [9j , it was technically chal- 
lenging to overlap PBG with ZnO gain spectrum. By 
adjusting the magnification of focused ion beam system, 
we were able to change the lattice constant a in fine steps 
of 3 nm over a wide range 100 — 160 nm. The ratio of the 
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air hole radius R to the lattice constant a was also varied 
from 0.20 to 0.30. In this way, we could tune PBG contin- 
uously through ZnO gain spectrum. We also introduced 
an artificial defect by missing an air hole. Structural 
analysis as in Fig. 1(b) gives the average displacement 
of a hole 5r ~ Q.22R. 

ZnO PhCS was optically pumped by the third harmon- 
ics of a pulsed Nd:YAG laser (A = 355 nm, 10 Hz repe- 
tition rate, 20 ps pulse width) at room temperature [g. 
In 8 x 8 /im patterns without intentionally-introduced 
structural defect, the ensemble-averaged lasing thresh- 
old exhibited a pronounced minimum at a — 115 — 130 
nm and R — 0.25a [Fig. 4(a)]. To understand this phe- 
nomenon, we calculated the photonic bands in ZnO PhCS 
using the computational technique described in Ref.[j|. 
The frequency dependence of ZnO refractive index was 
taken into account. In Fig. 4 (b), the wavelength of the 
dielectric bandedge Xd for the fundamental PBG of TE 
modes |§| is plotted against the lattice constant a. The 
structural parameters were extracted from the SEM of 
our samples. The ZnO slab thickness t = 180 nm, and 
R/a= 0.245. By comparing the lasing wavelength to A^ 
in Fig. 4(b), we confirmed that the lasing modes were lo- 
cated in the vicinity of the dielectric bandedge. This can 
be explained by two factors: (i) the electric field of the 
modes near the dielectric bandedge is concentrated in- 
side ZnO, and thus experience more gain [Fig. 3(c)]; (ii) 
the vanishing group velocity at the bandedge enhances 
light amplification The dip in the measured las- 

ing threshold [Fig. 4(a)] is attributed to spectral over- 
lap of the dielectric bandedge with ZnO gain spectrum. 
In Fig. 3(b), the measured lasing threshold decreases 
monotonously with the pattern size for a = 115 nm and 
118 nm. These data agree qualitatively with the numeri- 
cal simulation results shown in Fig. 3(a). In all patterns 
with intentionally missed air holes, the lasing modes were 
not pinned at the location of the missing hole due to the 
existence of better confined modes away from the defect. 
This observation is in line with our numerical simulation 
of large patterns with single artificial defect. 

In summary, the structural disorder may lead to self- 
optimization of optical confinement in a PhCS and forma- 
tion of high-Qt ot modes which serve as the lasing modes. 
In a sufficiently large PhCS with short-range disorder, 
a microcavity with balanced Q± and Q|| can be formed 
spontaneously without any carefully-designed structural 
defects. Despite the disorder, photonic bandedge effect 
enables us to efficiently extract optical gain and to fine- 
tune the lasing wavelength from 383 nm to 407 nm with 
sample-to-sample fluctuation of about 5 nm in ZnO PhCS 
lasers. 

This work was supported by the National Science 
Foundation under the grant no. ECS-0244457. 
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FIG. 4: (color online) (a) Experimentally measured incident 
pump pulse energy at the lasing threshold (averaged over five 
to ten samples) as a function of lattice constant a. (b) The 
circles are the measured mean wavelength of lasing modes, the 
error bar depicts the range of lasing wavelengths. The solid 
curve represents the wavelength of the calculated dielectric 
bandedge in ZnO PhCS as a function of a. The shade of 
the background qualitatively describes the position and width 
of ZnO gain spectrum. 
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